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Abstract 

Some old and new evidence for the existence of the string (planar random surfaces) 
representation of multicolour QCD are reviewed. They concern the random surface rep- 
resentation of the strong coupling expansion in lattice multicolour gauge theory in any 
dimension. 

Our old idea of modified strong coupling expansion in terms of planar random surfaces, 
valid for the physical weak coupling phase of the four-dimensional QCD, is explained in 
detail. Some checks of the validity of this expansion are proposed. (The lectures given in 
the Trieste Spring School and Workshop-1993 on String Theory) 



1. Introduction 



The principal model considered in this communication will be the Wilson lattice 
U(N) gauge theory without quarks, particularly its large N (multicolour) limit. 



The action of the theory is defined on the hypercubic D-dimensional lattice with 
the vertices labeled by the D-dimensional vectors x and the vectors of the links as 
(jl, u, etc. (so that the plaquette is denoted as (x, n,v))- 

S[U] = N(3tr £ U X:f ,U x+ ^U:+^U+„ + compl. conj. (1.1) 

X,fl,l/ 

The partition function Z, the free energy F and the Wilson average W(C) along 
a closed loop C on the lattice are now expressed as: 

Z = e N2p = f H ri [dU x ^} H expN(3S[U] (1.2) 

W(C) =< ^H MC U„, > (1.3) 

with the obvious definition of the average < ... >. 

The Haar measure [cIU}h for the U(N) integrals is characterised by the unitarity 
condition 

U+U = UU + = I (1.4) 

The famous paper of K.Wilson M, which explains the confinement of quarks in 
terms of the strong coupling (SC) expansion in lattice QCD, left a wide area in which 
to improve upon ; this lattice "superconfinement" resulted in a string tension K(/3) 
which was wrongly scaled to describe the weak coupling (WC) phase. Namely, from 
the Wilson calculations for the planar loop C of sufficiently large size one obtains in 
the (SC) limit the area law asymptotics: 

W(C) ~ exp[-K(f3)A min (C)\ (1.5) 
where A min (C) is the minimal area of the surface spent on C, with 

tf(/3)~-log03) (1.6) 

instead of the correct physical asymptotics 

K(/3)~A 2 exp[-^/3] (1.7) 

dictated by the dimensional transmutation mechanism of asymptotically free theory. 

It is widely believed, after years of computer simulations of the lattice SU(N) 
gauge theory, that in the case N=2 and N=3 gauge groups, if we move from small 



(3, where ( |1.6| ) is valid, to larger /?, we meet a sharp crossover: an abrupt change 
in behaviour in physical quantities as functions of /3, before we reach the physical 
confinement with the string tension ( |1.7| ). Nobody knows whether we can describe 
at least in principle the physical phase by summing up the SC series. Even if the 
answer is positive, practically, it turns out to be impossible. 

The situation becomes even worse if N increases. The crossover turns into a real 
first order phase transition, toppling any hope that the SC expansion can be extended 
into the physical phase. This phase transition corresponds to the breakdown of the 
Z(N) symmetry of the centre of the gauge group. It has indeed to be broken in the 
physical phase of the theory, since we know very well that each gauge matrix U on 
the lattice has to fluctuate close to some singled out value on the group, which leaves 
no space for the discrete Z(N) degrees of freedom. This phase transition is clearly 
seen even in rather naive mean field calculations. 

In the limit of large N, which is our only hope for any free string representation (if 
not exact solution) of QCD, a new complification arizes, namely, the Gross- Witten 
(GW) phase transition 0. This is also characterised by a narrow distribution of 
any particular lattice gauge variable around some singled out point on the group 
space, once one increases (3 beyond the GW transition point. If we consider the 
eigenvalues of a gauge variable (which are the gauge invariants) we can say that their 
distribution in the WC phase does not take the whole unit circle, as in the SC phase, 
but is squeezed into a smaller interval. The Z(N) symmetry of the group centre 
will be broken anyway. Note that although the neighbouring gauge variables should 
not be far away from each other on the group space, on the scale of the correlation 
radius of the theory (in other words, physical mass, or confinement scale), different 
U-variables cover the whole group space. This property itself can be viewed as a 
definition of the correlation radius. 

All this looks like very bad news for the SC expansion in the Wilson lattice theory. 
Even the hope of an exact string representation of QCD would be significantly reduced 
by this failure, since one of the few indications of its existence is the representation 
of the SC expansion in terms of the sum over non-interacting random surfaces on the 
hypercubic D-dimensional lattice, found in 0, [| and advanced in || ||, [F| (we will 
present below a short discription of this construction). 

Nevertheless, Wilson's original idea of confinement on the lattice seems to be too 
nice to discard so easily. A long time ago, the author proposed a possible way out 
of the impasse [|J. Namely, it was suggested to modify the SC expansion by the 
introduction of new weights on the links of random surfaces, which have a non-trivial 
dependence on (3 and must be calculated selfconsistently from the unitarity condition 
(|1.4j) for the gauge variables. These weights can be expressed in terms of averages of 



traces of products of the Lagrange multipliers for the constraints (|1.4j ) on each link. 
The weights do not fluctuate in the large N limit due to the factorization property, 
and are uniform in the physical space in light of the translational and rotational in- 
variance. So, they can be considered as the effective "world sheet coupling constants" , 



though this notion is a bit vague on the lattice. 



The whole construction is very close to the Stanley solution of the N-vector field 
theory (on the lattice) in the large N limit (though our construction is as yet far from 
a solution). Let us recall these simple and nice arguments of Stanley, to compare 
with ours. 

The partition function of the N-vector field n a (x) normalized as 

nl(x) = 1 (1.8) 

is defined as follows: 



Z = I H x D N n(?)8(?4(x) - 1) expLY^ £ n b {x)n b {x + //)] (1.9) 



We can try to investigate this model in terms of the SC expansion by expanding ( |1.9j ) 
in powers of (3 and integrating order by order with respect to the compact measure in 
n(x). The strong coupling diagrams will look like sums over paths with the weights 
^length ^ w j 1 j c j 1 m ight have nontrivial couplings in the points of selfintersections. Let 
us demonstrate, say, the picture corresponding to the constraint (|1.8| ): 



<K(x)>=i (i.io) 

The typical strong coupling diagram is shown in fig.l. In the large N limit the 
diagrams will simply be trees of loops attached to each other at single points in 
physical space. 

Let us remind ourselves how to sum up these trees. 

We introduce the lagrange multiplier field a(x) to impose the constraint ( |1.4| ) on 
the integration measure, and then integrate out the n-field to obtain the following 
functional integral: 

f N 
Z = J U x Da ( x ) e MN^a(x)--Tr\og(~(3(6 XtX+l , + 5 x , x ^ + a(x))} (1.11) 

In the limit of large N we can solve this problem by looking for the translationary 
invariant saddle point for a, 

^# = (1.12) 

6a(x) 



which gives the equation: 



2tt I 

, . d 2 p— — (1.13) 

In the limit (3 — > oo we shift the variable a 

a = [3(4-m 2 ) (1.14) 
expand the cosj? in (|1.13|) (since we are going to get the continuous limit) and obtain 

(zttj 2 J p 2 + m 2 47r m 2 
which gives the well-known dimensional transmutation for the physical mass scale: 

where A is the ultraviolet cutoff. 

The end result is a multiplet of N scalar noninteracting particles with the same 
mass ( [L.16 ). This result is exact in the large N limit. One can easily see that the 



saddle point condition ( |1.13|) corresponds to the normalisation condition ( 1.10Q . 
If we expand the formula ( |1.13| ) in powers of /3 we get: 

J _ ^ @ ^Length{P xx ) (1 17) 

pathsP xx ^ a 

where P xx is a path starting and finishing at the same point x on the lattice. Instead 
of summing over the whole tree of the loops we sum over only one free closed path 
with the renormalised weight: 

/ ft \Length(P xx ) _ _ ^\Length{P xx ) q 

< a > 



We see that this sum over random paths diverges at m = 0, as it should do, in 
order to have an appropriate continuous limit. 

This is an example of a modified SC (MSC) expansion which works in the physical 
SC phase. 



It will be useful for the future comparison with QCD to show the relationship be- 
tween the lagrange multiplier field and a vector field condensate (I thank A.Polyakov 
for his comments on this subject). We have the following equation for the correlation 
function of the vector field in the background of the lagrange multiplier: 

(3 < n b( x + v) n b(y) > + < an b (x)n b (y) > +8 Xty = (1.19) 

/x=±l,±2 

For x = y we obtain in the continuous limit: 

< a >= Af3 - pA- 2 n 2 < n b (x)An b (x) > -1 (1.20) 
= 4/3 + ^/^-1 (1.21) 
= 4/3-f^log^ (1.22) 

This is the relation between the modified weight in the sum over paths and the vector 
field condensate. Of course, it agrees with ( |1.16|) . 

We presented here this well-known solution of Stenley model since it will be our 
guideline for the construction of the MSC expansion in lattice QCD, which looks to 
be generalizable to the physical WC phase. 



2. Free random surface representation of the strong coupling 
expansion in D-dimensional lattice gauge theory 

In this section we are going to reformulate the standard SC expansion of Wilson 
gauge theory in terms of FREE random surfaces on the D-dimensional hypercubic 
lattice. This formal representation will be exact order by order in (3. The result will 
consist of a description of the elementary geometrical objects (plaquettes, "saddles" 
et.c) from which we construct this surface. 

To demonstrate the idea let us start from a simpler model of random surfaces on 
the lattice: the Weingarten model. It has the same action as Wilson gauge theory, 
but instead of the unitary Haar measure one takes a gaussian measure: 

[dU] w = d 2N2 Uexp[-NtrU + U] (2.1) 

It can be represented in terms of lattice random surfaces by the standard SC expan- 
sion. We expand the exponent in ( |1.2|) in powers of (3 and perform gaussian integrals 
using Wick's theorem: 

J[d 2N2 U} w U±U kl = ^5 u 5 3k (2.2) 



r 2 1 1 

J [d 2N U} w U+ ji U klh U + i 2 j2Uk 2 i 2 = j^S ilh S jlkl 5i 2h S j2k2 + j^S h i 2 5j lk2 8 i2k 8 j2kl (2.3) 
and so on. 

Geometrically this means that we can glue plaquettes together by means of these 
integrations over the common link variables, as is shown in fig. 2 for (|2.2|) for two 
plaquettes, and in fig. 3 for (|2.3|) for four plaquettes. Continuing this process we 
arrive at closed planar surfaces built on the hypercubic lattice. Due to the matrix 
structure of the theory, every connected piece of the surface will be weighed with the 
weight N 2 ~ 2G , where G is its genus, this factor is explained in the standard way, for 
matrix models. In the Weingarten model the 'tHooft limit of large N leads precisely 
to planar surfaces. Due to the equal weights of different terms in the formulae (|2.2j )- 
( |2.3| ), there are no extra weights at the intersections of surfaces which means the 
surfaces are the world sheets of a "free string" . 

Finally, the SC expansion for the free energy in the Weingarten model can be 
represented in the large N limit as a sum over planar surfaces aw on the hypercubic 
lattice: 

F(/3) = ^(3 A M (2.4) 

aw 

where A(aw) is the area of the surface (the number of plaquettes from which it is 
built). 

It is known that for any dimension higher than 1 this model has a pathological 
behaviour and degenerates into tree like configurations at the critical point M. 

What will change if we take instead of (|2.1|) the U(N) Haar measure of the Wilson 
gauge theory? Nothing will be different for the gluing of two plaquettes, since the 
integral ( p.2|) will be the same. However, for four plaquettes, we will have instead of 

(U): 

I[d 2N } H Ut lh U klh U+i 2 j 2 U k2l2 

— jV(AT 2 _i) [^hh^kih^hh^k^ji ~~ ^hl2^k2j2^i2h^kiji] (2-6) 
= Jp8i 1 i 1 5j lkl $i 2 l 2 5j 2 k 2 ~ Jp^hh^k 1 j 2 8i 2 l 2 S k2 j 1 (2-7) 

+ jpifiiihfijikifii 2 i 2 fij 2 k 2 — ]^5^ih^fcij2^2/2^2ji + ••• + (h, ji^h, 32) (2.8) 



The first line in the r.h.s. of (2^8) represents the original Wick contractions of 
whereas the second corresponds to new, cyclic contractions of indices (and plaque- 
ttes), etc. 

In the large N limit, we have to expand the N-dependent coefficients in (|2.8| ) in 
powers of 1/N and interprete each term geometrically in terms of pieces of a free 
planar random surface on the D-dimensional hypercubic lattice. The most natural 



interpretation is presented graphically in fig.4. The first term represents one of Wick 
contractions of indices (the rest of the gauge variables on these four plaquettes is 
contracted into two separate traces at the boundaries of two disconnected couples of 
plaquettes) . 

The second term corresponds to the cyclic coupling of indices (the rest of the 
gauge variables are contracted in a single trace around a single boundary). This 
term looks like a saddle (though quite a singular one, with zero radius of curvature). 
Topologically it is a disc which is glued into the random surface. It has the an 
extra 1/N power with respect to the previous term, corresponding to smaller Euler 
characteristics k (it is equal to 1 for each disc). In this way the topological expansion 
of 'tHooft attaches to our surfaces the factor: 

N 2-2G = n k 

The saddles seem to play an important role in the whole construction. Note that 
the sign is negative in front of this term. 

The last term in fig. 3 represent the next term of the 1/N expansion of the coeffi- 
cient in the first term in the l.h.s. of fl2.8| ) and can be described as a tube connecting 
two couples of plaquettes. It again has the topology of a disc. 

All these terms can in principle contribute in the large N limit. Higher order 
terms in the 1/N expansion in ( |2.8|) describe higher order topologies. 

This procedure can be continued for higher n-correlators of gauge matrices of the 
type 

/ [d 2N2 \ H Ut lh U kxh Un nJn U knln (2.10) 

These correlators, being calculated in the same manner, give all possible connected or 
disconnected objects like multiple "saddles", drawn in fig. 5, and give rise to a cyclic 
contraction of indices in ( [2.10|) , like 

8i 1 i 1 5k 1 j 2 5i 3 i 3 ....8 kn j 1 (2-H) 

as well as tubes and their mixtures, which can be build from 2n plaquettes. Every 
connected part of such a configuration is accompanied by some numerical coefficient 
which corresponds to appropriate index contractions in the integral ( |2.10| ). 

This interpretation of the strong coupling expansion in terms of random lattice 
surfaces was proposed in || |]. The factors corresponding to the n-saddles, consisting 
from cycling gluings of 2n plaquettes, were found there to be the Catalan numbers 
f n (see next section for their calculation): 



(2.12) 



so that 

A = l, / 2 = -l, ,/ 3 = 2, / 4 = -5, / 5 = 14 (2.13) 



The whole variety of coefficients corresponding to this zoo of objects was calcu- 
lated in ||. Let us comment also that the rules for the sum over surfaces found here 
for the free energy, are directly generalisable to the Wilson average W(C). One can 
consider the gauge variables forming the loop factor in (|1.3|) as the edge of a surface, 
to which the plaquettes can be attached either by Weingarten type contractions, or by 
saddles, tubes, et.c. The Wilson average can be viewed as an open string amplitude. 

Fortunately, as was shown by Kostov || [7], with the exception of the multiple 
saddles, these complicated objects can be ignored and the sum over surfaces can be 
reduced to a sum over one-link-irreducible surfaces with multiple saddles. It was 
shown by Kostov that any surfaces which can be cut into two pieces by cutting along 
a single link (one-link-reducible surfaces) cancel each other due to sign changing 
factors. 

This theorem of Kostov is a direct consequence of the XSunitarity of the gauge 
variables ( |1.4|) . Let us sketch the proof of it (the details can be found in (|7j]). 

One can use the "backtracking" condition for an arbitrary Wilson average ( |1.3| ): 
if we cut the contour C in a point and glue in this cut a path consisting from two 
links I * l^ 1 going there and back in the same direction \x on the lattice, the Wilson 
average will not change (see fig. 5): 

W(C*l*r l ) = W{C) (2.14) 

Now consider dressed saddles: by definition these correspond to the sum of all surfaces 
attached to the link on which the saddle is situated. Namely, it is the bare saddle, 
considered above, plus all surfaces connected to it by "tubes" (also considered above). 
We can have in the sum over surfaces in the l.h.s. of (|2.14j) three situations: a dressed 



link vertex of the n-th order can be attached to both links I, l" 1 , or two link vertices 
of the orders (n-k) and k can be attached separately to each link, or there will be 
no plaquettes attached to these links at all (which corresponds precisely to W(C) 
without back tracking). So, introducing the weights Fk of dressed saddles, we obtain 
from (glil) : 

n 

J2 W n {C)[F n + Fn-kF k } + W(C) = W{C) (2.15) 



fc=0 



where W n (C) is the sum over surfaces spanned on the contour C with dressed n-saddle 
attached to two backtracking links. 

Since C can be arbitrary here, we conclude that 



n 

F n + y"l F n _kFk = 

fc=0 



(2.16) 



which defines the Catalan numbers ( [2.12 ). From here we conclude that the weights 



of the dressed saddles are also equal to the Catalan numbers: F n = f n . Hence, 
we can throw away the one-link- irreducible surface, together with the objects more 
complicated then multiple saddles, from the sum over surfaces. 

As usual, the partition function corresponds to the sum of various disconnected 
random surfaces, each having the topology of a sphere. Since these surfaces are non- 
selfinteracting (there is no excluded volume problem for them) we can be sure that 
the free energy corresponds to the sum over only connected random surfaces. 

Finally, the random surface representation of the SC expansion in the Wilson 
gauge theory looks quite simple: 

F(f3)=J2P Areai(T) Ilfns (2.17) 

cr sa 

where a are the one-link-irreducible planar surfaces built from plaquettes glued to- 
gether in saddles denoted by s, on the D-dimensional hypercubic lattice. f ns are 
the Catalan factors attached to these saddles (we consider normal contraction of two 
plaquettes saddle of order 

The analogous representation for W(C) can be given in terms of a sum over 
surfaces spanned on a loop C. Very few things are known about the critical behaviour 
of this SC expansion in multicolour Wilson gauge theory. The no-go theorem of || is 
not directly applicable here, since we have sign changing terms. It might or might not 
have a nontrivial critical behaviour, but its direct continuation to the physical WC 
phase seems quite unprobable in virtue of one of the two phase transitions mentioned 
in the introduction. 

Can we modify the random surface representation ( |2.17| ) in order to describe the 
WC phase of QCD as well? We will propose a possible modification in the next 
section, and then we will give some arguments in favour of it. 



3. Modified strong coupling expansion for the physical phase 

in QCD A 

To explain the idea of the modified SC expansion, proposed in || and tested 
on some examples in ||, we consider two similar, but technically slightly different 
approaches. 

The first one, originally proposed in j3| , leads to a simpler geometrical picture for 
random surfaces, whereas the second based on the representation introduced in |7| 
leads to a more promising quantitative scheme. 



3.1 First construction for modified SC expansion 



In the first case we introduce for each link gauge variable a hermitean lagrange 
multiplier matrix a M , parametrizing the Haar measure in one of two following ways: 



It is our choice whether to take the second or the third line as a definition of a for 
any link variable. Let us choose it in the most symmetric way: we classify all the 
vertices of the D dimensional lattice in chess order, as even and odd, and we take the 
definition of according to the second line of ( |3.3| ), if the corresponding link goes 
from an even to an odd vertex (in the positive direction of the coordinate axes fi), 
and the line three of (13.30, if otherwise. 



Now we have a double functional integral for (O): in hermitean matrices a and 



in complex matrices U. Let us perform first the integral over Ws by means of strong 
coupling expansion in f3. In each order the integrals are now purely gaussian, with 
the matrix propagators equal to 



/ 



d™ 2 e m ™ u+u U+U kl = ^(a-%6 jk (3.4) 



for the "even" links, and 



[dU] H = d 2N2 U5 N \U + U -I) = 
fd N2 ae m ™V- uu+ ) 



(3.1) 
(3.2) 
(3.3) 



/ d™ 2 e m ™ uu+ Ut 3 U kl = ^5 a {a-% 3 (3.5) 

for the "odd" links. 

The integration over Ws can be performed according to the same rules as for the 
standard Weingarten model with the measure fl2.ip . We will get the sum over exacly 
the same hypercubic surfaces aw as in Waingarten model, but since the propagators 
are now modified, new a-depending weights should be attached to every even vertex 
es of the surface: 

F(p) = j2P Area(CTw) n mM ( 3 - 6 ) 

uw esyya 

where 

SV} = (0) =< j^ tl [ a ^ a ^-- a ^Y 1 > ( 3 - 7 ) 

where fix, //2, £*n are the directions, either positiv or negativ, of links around a 
given even vertex on the surface). 



Instead of futher integration over a's, we have already substituted in ( |3.7| ) the 
emerging traces of products of a;~ 1 's by their averages, according to the large N fac- 
torization theorem [II], [□]]. We see from this construction, that we can formally 
represent any physical quantity in multicolour lattice QCD as a sum over random 
hypercubic surfaces (Weingarten type surfaces in this case) with special factors at- 
tached to the vertices of surfaces. In virtue of translational and rotational invariances 
of any physical averages, the (/-factors should not depend on the position of the ver- 
tex of the surface in the D-dimensional lattice, but only on the sequence of the links 
fix, around a given vertex of the surface (up to the obvious rotations). So 

play the role of " string coupling constants" , though this notion is quite vague for the 
lattice strings. 

These constants are in fact non-trivial functions of (3 and should be calculated 
separately. One can propose the following formal method for it. 

In the large N limit, we can in principle define all g-f actors from the chain of 
obvious equations, following from the unitarity of gauge matrices: 

< ^K U ^- U t U ^ >= 1 (3-8) 

where //i, ...,//„ is any sequence of directions of links around an even vertex on the 
lattice. The l.h.s. of ( |3.8|) is just a wilson average for a loop C{n} consisting of 
only backtrackings along the links surrounding a given vertex. If we re-express these 
Wilson averages in terms of the same sums over random surfaces spanned on these 
contours, as it was done for the free energy (|3.6|), we obtain the following conditions 
on g's: 

E P Area((7w) U = l (3.9) 

a w ,da w =C{^} es °W 

where the sum is taken over all the surfaces having as a boundary the abovemantioned 
backtracking contour C{ji}. 

This chain of nonlinear equations can serve at least in principle for the calculation 
of the string couplings g^{(3). Note the similarity (not accidental) with the equation 
(|1.17|) for the N- vector field: there the sum over paths pinned to a point in the physical 
space, with the renormalized hopping parameter (3/ a, was equal to 1, as a consequence 
of the normalization of the vector field, where as here the sum over random surfaces 
pinned to a sequence of links around a vertex on the lattice, with renormalized string 
couplings, is equal to 1, as the consequence of unitarity. 

Of course, technically all this looks too complicated. There are too many factors 
g to calculate. Therefore we use this approach only to demonstrate that the sum over 
random surfaces is quite a natural representation of the Wilson multicolour QCD. 

Why should this modified strong coupling expansion work in pnysical WC phase 
of the theory? Our hope in this approach is that the equations (|3.9|) can have two 



different branches of solution: the SC branch, which corresponds to the standard (3 
expansion considered in the previous section and valid up to some critical j3 c , and 
the physical WC branch valid beyond f3 c . It is obvious that in the SC phase the g- 
factors will be just /^-independent numbers, where as in the WC phase they should be 
nontrivial functions of (3. The traces of a-matrices serve here as an order parameter 
for the corresponding phase transition (of Gross- Witten or Z N breaking type). We 
will show in the next section that this modified SC expansion works indeed in both 
phases in the simplest example: the one plaquette model. 

This construction is hardly useful for practical calculations, but the existence of 
a simple random surface representation of Weingarten type suggests that the search 
for the continuous QCD-string is not in vain. 



3.2 Second construction for modified SC expansion 

Let us describe the second construction for the modified SC expansion. The corre- 
sponding surfaces in this construction will be more complicated, including now the 
n-saddles. However technically it will be much more tractable since we will effectively 
obtain a sequence of weights g n (3 (traces of only one a-matrix), attached to n-saddles 
labeled by only one integer n. 

According to the trick proposed in (@), let us double every gauge variable on 
each link I: 

Ui -> XJiVf (3.10) 

We can integrate now in each of the matrices separately. This is the same Wilson 
gauge theory, since the matrix V can be easily absorbed into U, in virtue of the 
invariance of the Haar measure. 

Now the Haar measure on every link can be parametrized as: 



[dU] H [dV] H = J d^a.J d N2 a 2 Jd 2N2 U J d ^ Ve mra l{ I-U + U) e Ntra 2{ I-V + V) 

(3-11) 

We are ready now to integrate out the U, V variables by means of the formal (3- 
expansion with fixed a's. Since each of the link variables is doubled we have to apply 
the Wick theorem separately to U-half-link and V-half-link independently. Let us 
consider two examples of one link integrals over U and V with a fixed: with two 
link-variables: 

< (UV + ) ij (VU+) Kl >J / 0) y = ^6 a 8 jk jja- a (3.12) 
and with four link-variables: 

< (VU + ) lin (UV + ) klh (VU + )i 2 j 2 (UV + ) k2h >® v (3.13) 

= a?2 (%0£~ 2 ) 2 8i lh 8j lkl 8 i2 i 2 5j 2k2 (3-14) 



+ Wjj a 4 $hl2$k 2 j2$i2h$k 1 j 1 (3.15) 

+ 31^2, 32) (3.16) 
where we introduced the hermitean matrix 

1 1 

a 2 = Oi1oL20i\ (3-17) 

So the cyclic contraction of indices corresponding to a saddle, appears here in a 
natural manner, even on the phase of gaussian integrations. To every n-saddle a 
factor 

9n(P) =< ^a- 2n > (3.18) 

should be attached. As in the previous construction, in virtue of translational and 
rotational invariance, and the large N factorization theorem, we can already take the 
average for each factor g n , which will not depend on the position or orientation of 
the link. So g n 's are nontrivial functions of j3, labeled by only one integer n. The 
corresponding sum over surfaces a* will consist from surfaces built from plaquettes 
glued together by saddle-like configurations s: 

w)=E^ rea(CT * ) n^ (3-i9) 



Note that unlike the surfaces a in ( 2.1 7| ) the surfaces a* emerging here can be one 
link reducible. 

The random surface picture here is more complicated then in the previous con- 
struction, and a good question is whether we can describe the emerging saddles as 
special vertex operators in a continuous string theory. We will discuss this possibility 
in the concluding section. 

Now we have to calculate the weights g n . The best way to do it is to write down 
the effective action for two matrix variables 0.1,0.2 (as in the previous construction, 
they will not depend either on the space coordinates, or on orientations in the large 
N limit) defining these factors through eq. Q3.18| ): 



S e ff(ai, a 2 ) = iVtrfai + a 2 — log[aia 2 ]] + N 2 F[a l a>2\ (3.20) 

where F^ia^] is the same as F in ( |3.19| ), but with unaveraged traces of a's instead 
of # n 's. 

The partition function of QCD can be written as a two matrix problem: 

Z(P) = J d N2 ai J d N2 a 2 e s ^ ai ^ (3.21) 



We propose here to reduce the double integral in a±, «2 to an integral over the only 
matrix variable ( |3.17| ). For this purpose one has to integrate out the extra degrees 
of freedom. Introducing the variable (|3.17|) and 



7 = "i 



(3.22) 



we rewrite (|3.21 ) as 



Z(J3)= f d N2 a 2 J rf 7v2 7 exp(Artr[7 + 7 - 1 a 2 -log7] +N 2 F[a 2 }) (3.23) 



Let us compare it with the well-known integral over the unitary matrix U in an 
l 

external field a 2 : 



N 2 B[a] 



By introducing the parametrization 



[dU] H e 



Ntra 2 [U+U+] 



[du] H = d 2N2 u5 N \u + u -i) = J ^V 7 ^ 7 -^ 



(3.24) 



(3.25) 



and integrating out U we can see that the integral over 7 in ( [3.231 ) is precisely given 
by (PI) : 

Z(J3) = J d N 'a 2 exp N 2 (B[a] + F[a 2 }) (3.26) 



One can find in [|1^] the expression for both the SC and WC phases of this integral. 



We need only the WC branch of this solution, found in fll"3||: 



N 



1 3 

B[a] = 2/N «fe - 2^72 12 l °9l a k + a j) ~ 7 

k=l k>j 



(3.27) 



where ai, ...ajv are the eigenvalues of the a-matrix. 

One can easily obtain this result by using the representation U = e tA — 1 + 
%A — I A 2 + ... of the unitary matrix in ( |3.24|) and keeping only guassian terms. This 
approximate calculation gives nevertheless the exact result in the WC phase, as is 
the case for the matrix integral in 



Finally, we can define the momenta ( |3.7| ) as the saddle point condition on the 
density of eigenvalues p{a) in the one- matrix integral ( |3.26|) : 



00 „ 1 1 

1 + £ a - 2n ~ l F n [ gil g 2) ...] = P / dpp(p)( + i— 

~[ J v « — p a 4- 



(3.28) 



where F n are the sums over the surfaces built according the same rules as ( |3.19| ), 
but attached to a fixed saddle of the order n. So to compute g n 's we have first to 
calculate these sums over surfaces (the most nontrivial part of the problem) then find 
p(a) from ( |3.28p as a function of g> n 's and then to find them from the selfconsistency 



condition: 

P dap(a)a~ 2n = g n (p) (3.29) 



One can recognize the geometrical similarity of the eq. ( |3.28| ) with the eqs. ( |1.17| ) 



The idea of the modified SC here is similar to the previous construction: we hope 
that the factors g n {(3) will have two different branches as functions of (3. In the SC 
phase they are just /3-independent numbers, where as in the physical WC phase their 
behaviour changes and they turn out to be nontrivial functions of (3. We will find 
these two branches for the simple one-plaquette model in the next section, and then 
we will demonstrate their existence in the 4D Wilson theory. 



4. Simple examples: one-plaquette model 

In this section we will check the idea of the modified SC expansion on a simple 
example: one-plaquette model, which was defined and solved in 0. Its partition 
function is: 

Z P (0)= J[dU] H e Nt ^ u+u+ ^ (4.1) 

We will use a method for its solution explaned in II which recalls our first construc- 
tion. Namely we parametrize the Haar measure in ( |4.1| ) by the lagrange multiplier 
a, as in ( |3.3p , and after the integration over the complex matrix U we obtain the 
following 1-matrix model: 

Z P (P) = [ d N2 ae N ^ a - losa+l32a ~^ (4.2) 



By standard methods |L5[ we obtain the integral equation for the density of eigen- 
values p(a): 

-l + - + P 2 a~ 2 = P fdfipUt)—^— (4.3) 
a J a — p 

If we introduce the analytic (outside the cut) function 



F(a) = P f dpp(p)^- (4.4) 
J a — p 



which satisfies the condition: 

F(a) -^oo 2/a (4.5) 

and has no singularities at a = 0, with its imaginary part being equal to irp(a), 
we obtain two solutions for it, separated by the Gross- Witten phase transition at 
P c = 1/2. The SC solution is valid only for f3 < 1/2 and corresponds to the density 

p( a ) = ^(l + I) (4.6) 

The cut is collapsed to a point here, and the calculation of momenta of a should be 
understood as a contour integral around the origin. So we obtain: 

<^a~ n >= j dap(a)a- n =1, if n = 0, 1 (4.7) 

= 0, for other n (4.8) 

It easy to obtain the following expression for the plaquette averages: 

W n ((3) =< ^(U n + U +n ) >= 2/T < ^cT n > (4.9) 

We conclude from the two last formulae that: 

W 1 ((3) = (3 (4.10) 

and all other W n are equal to zero in the SC phase, as it should be. 

In the WC phase {(5 > 1/2) we find the solution for F(a) in the form 



F(a) = -1 + _ + £- + (_ + 4) V / -(« 2 + 2(l-/5)a + /3 2 ) (4.11) 
Comparing it with the expansion 



oo t 

F(a) = £ a k < -a-"- 1 > (4.12) 



we can obtain: 



W 1 = 1 - h (4.13) 



4/3 

2^ 



^ 2 = (l-^) 2 ,etc. (4.14) 



We see that unlike the SC phase, all W n are functions of ft tending to 1 as ft — > oo, 
as it should be. 

It is easy to see that in the SC phase the large ft asymptotics of the momenta of 
a are always: 

< -^a- n >^oo/T n (4.15) 

which is similar to the asymptotics of a single lagrange multiplier of the vector field, 
as one sees from Q1.14 ). 



It is curious to see how the unitarity relations (U + U) n = I are satisfied here. Say, 
from the integral over U with fixed a we have 

< ^(U + U) >=< > +ft 2 < |a" 2 >= 1 (4.16) 



One can observe that the WC solution and the SC solution satisfy (|4.16 ) in very 
different ways. 

Of course the example considered here is the simplest possible one: the modified 
SC expansion consists here only from one term: in ( [Op this is the last term in the 
exponent, where the SC expansion parameter ft 2 is modified by the weight < ^a" 1 >• 
In the SC phase this weight is just one, and ft 2 corresponds to the area of the "surface" 
consisting from two plaquettes. In the WC phase the weight is already nontrivial, as 
you see from (|4.14|) . The same is true for the "Wilson loops" W n (ft). 

One might think that the situation considered in this example is too simple to 
test the validity of the modification of SC, since we have only one surface in the sum 
over lattice surfaces here. Let us try another formulation of this model which already 
has an infinite sum over the " surfaces" (built on one single plaqutte from an arbitrary 
number of copies of this plaquette). Namely let us test our second construction of 
the previous section. The doubling of the matrice U gives instead of (O) the two 
unitary matrix integrals: 

Z P (ft) = J [dU] H J [dV} H e m ^ uv++vu+ ^ (4.17) 

Introducing the lagrange multipliers as in (|3.11|) and integrating over U, V we obtain 



the sum over surfaces with modified weights in (|3.26 ) in the form: 



F(a) = -| log(a 2 - ft 2 ) = - log a 2 + £ ^|«" 2fc (4.18) 



The saddle point equation (|3.28 ) now reads: 



- 1 + 



a 



a 2 - ft 2 



P f dnp{n)(-^— + \-^- 
J v ct — (j, A a + /i 



(4.19) 



This equation seems to be rather nontrivial, but nevertheless it can be easily 
solved in the relevant WC phase. Let us shift the matrix variable a, defining a new 
variable e: 

a = (3 + e (4.20) 

The eq. ( |4. 19[ ) now reads: 



e + (3 



p /^b + hsb) (4 - 21) 



2ef3 + e 2 J ry ^e-e' 2 2f3 + e 

We can try now to solve this equation by expanding the second term in the l.h.s. 
and the second term in the r.h.s. in powers of /3~ x . In the r.h.s. we will obtain 
in this way the momenta < e > as the coefficients of expansion. It happensto be 
self consistent to set all of them to zero. 

<e n >=0, for n = 1,2,3,... (4.22) 

Then the equation ( (4.21| ) turns out to be 

-l + L = P J de ' p{e ' )7 L_ (4.23) 

with the solution: 



As in the eq.( [4.6|) the corresponding integrals should be understood as contour inte- 
grals around the origin. 

It seems to be strange that the solution ( [4.24D does not depend on (3 and is very 



similar to the SC one (refscr). Nevertheless it reproduces the correct WC behaviour 
of Wilson averages, say, 

Wl <fi> =< ±<uv+ + w+) >=< >= f ieAe) _L_ = 1 -.L (4 . 25) 

We see from the last formula, that the modified SC expansion contains here an 
infinite number of terms: 

Wl ((3)=(3j2(3 2 k<-a^ k - 2 > (4.26) 

k=0 



At (3 = oo this series diverge since a ~ j3. This divergence is situated at e = 
0, therefore it gives a necessary contribution to the integral in ( |4.25| ). For every 



particular term in the expansion ( |4.26| ) we would loose this singularity, therefore 
to get the correct weights we have to regularise the singularity in the density of 
eigenvalues ( |4.24| ). 

One of the regularizations would be the above mentioned expansion of ( |4.21| ) in 
1/(3 up to a finite number M of terms and keeping the corresponding momenta ( |4.22| ) 
nonzero for a while. Then we can solve the resulting equation by standard methods 
and define the inverse momenta < jjOt~ k > (the weights of the "surface") from the 
self-consistency condition. The " sum over surfaces" with these weights will reproduce 
the correct result ( }4.25|) for sufficiently big M. This method was checked numerically 
in || on a similar construction and the convergency with M turned out to be very 
fast. 

Anyway, this example shows that the modified SC expansion works well in the 
WC phase of the one-plaquette model and obeys the following general features which 
should hold in more realistic models: 

1. There exist two different branches (SC and WC) of the solution for the weights 
of the modified SC expansion: they are numbers in the SC phase, and non-trivial 
functions of (3 in the WC phase. 

2. The matrix lagrange multiplier has the asymptotic behaviour: 

a ~ 0, for -> oo (4.27) 

tending to cancel the big factor p Area in the sum over surfaces. 

3. The corresponding sum over surfaces diverges at the point P = oo and needs to 
be regularized. 



5. The existence for the modified weights of random surface in 

physical phase of QCD4 

The problem of calculating the sum over surfaces in a physical four dimensional 
situation in the formulae ( |3.19| ) is non-trivial. Even its reformulation in terms of a 
continuous functional integral of some string theory with a definite two dimensional 
world sheet action is far from being achieved. 

Our task here is rather to demonstrate the possibility of the existence of such a 
representation in principle, then to get a well elaborated quantitative approach to 
the multicolour QCD. 

Namely, we will show the relationship between gluon condensates and the mo- 
menta of the lagrange multiplier matrix. 



Let us vary the action ( |1.1Q with the Haar measure parametrized as in ( |3.11|) with 
respect to the fields U + 5U + , V5V. We get the following matrix operatorial equations 
of motion: 

= PJ2( uv+ )*Auv + ) x+ ^(vu + ) x+ ^(vu + ) x ^ (5.1) 

V 

and 

a 2 x lAt = PY,( UV+ )*AUV + ) X+ ^(VU + ) X+ ^(VU + ) X , U (5.2) 

V 

Since we know that any physical gauge invariant quantity depends only on the com- 
bination a.ia.2 of these matrices, we find from last equations: 

a = ol\Ol<i = (5-3) 
(3J2 UUU2 (UV+) X , U1 (UV+) X+ ^(VU+) X ^ 

(5.5) 



The product of matrices in the r.h.s. runs along the 6 link boundary of two 
plaquettes attached to each other along a fixed link //. This "chair-like" configuration 
of plaquettes was used extensively for various improvements of the Wilson action in 
lattice computer simulations of QCD. 

So any momentum of the matrix a which acquires the space and direction inde- 
pendent average in the large N limit, can be represented as an average of the power 
of the combination in the r.h.s. of (|5.5f) . Any positive power looks like a saddle 
configuration of the n-th order (built from 2n plaquettes) presented in fig. 4: 

<^a">= !^<%(EwUwU a ^JJ^U^wU^JJ^) n > (5.6) 

+ contact terms (5.7) 



If we represented the n'th power in the r.h.s. as a sum of individual terms it 
would be expressed through Wilson loops running around the boundaries of saddles 
of the n-th order (we restored the standard Wilson gauge variables U, without the 
doubling on each link). 

For negative n's ( |5.7| ) gives a local expression of the weights of the modified SC 
expansion in ( |3.19| ). Again, the situation is very similar here to the calculation of 
the lagrange multiplier through the local condensate of the vector field, give in the 
Introduction. 

Let us discuss a possible continuous limit for the expression ( |5.7|) . it is clear that 
for (3 — > oo any wilson loop tends to 1, so, a — > 4(D — l) 2 /3. As usually, we make a 
shift: 

a = A(D - if (3 + e (5.8) 



It already shows that the big weight f3 Area will be cancelled in modified SC expansion 
( |8.19| ). Going to the local limit in ( |5.7| ) we obtain 

< a_ n > (5.9) 
~ ([2D - 2}-^l - < e > +...) (5.10) 

~ (3 < jjF^ v (x) >qcd> + substractions (5.11) 

(5.12) 

where F^ u = A u , IM A fl , u + [^4^, A u ] is the gauge field strength. 

We see that the weights of the modified SC expansion are defined through the 
local condensates of the gluon field A^(x) in QCD and can in principle be calculated 
perturbatively, from the Feynman perturbation technique. All terms proportional to 
the powers of cut-off in ( |5.12|) should be substracted, since they should cancel with 



the contact terms, as happend for the vector-field. 

However, the perturbative calculation cannot give us the most interesting expo- 
nentially small terms in the gluon condensates, which should in principle define the 
renormalized string tension. This problem cannot be solved if we are unable to cal- 
culate the corresponding sum over surfaces ( |3.19|) and get the effective action for 



the a variables. The only purpose of the relations found here, between the gluon 
condensates and the momenta of the lagrange multiplier, was to show the existence 
of the new nontrivial weights in the modified SC expansion leading to a lattice string 
representation for the multicolour gauge theory. 



6. Conclusions 



Presenting this lattice construction of the random surface representation in QCD, 
we did not want to create an impression that a real chromodynamical string theory 
is around the corner. All of this can serve only as an intuitive picture from which 
one can probe different continuum string world sheet lagrangians. We tried here 
to demonstrate in principle the existence of the modified strong coupling expansion 
with the weights of the expansion being defined selfconsistently in the physical weak 
coupling phase. We managed to demonstrate the existence of these weights in a 
toy one-plaquette model, and then we gave some hints as to how to calculate them 
through gluon condensates. 

What are the main ingredients of this world sheet lagrangean? 

As it was noted in papers and pOfl , the fact that the Wilson loop average 
should obey the Makeenko-Migdal equations implies the existence (finiteness) of the 



area variational derivative of the Wilson loop |19| . This means that the cosmological 
constant term should be absent, and the action contains the term: 



/ d^T^daX^X, (6.1) 

where A M (£) is the world sheet coordinate in the physical space, and is either the 
independent antisymmetric field on the world sheet (in [TSJ), or an antisymmetric 
function of X M (in g). 

This might be a good starting point but it is far from the final action. According 
to our construction of this paper, a probable object to be added with a new coupling 
constant to this action would be a saddle like configuration. This saddle is quite 
singular: for example, for the saddle of the second order shown in fig. 3, we have a 
deficit of the angle (curvature) equal to Am (since it consists from 4 plaquettes) and 
concentrated in one point on the on the world sheet. The radii of curvature at these 
points are equal to zero. Nevertheless, this singular object is localized on the world 
sheet and can in principle be described by a local operator in the action, depending 
on the extrinsic geometry. The problem is to guess the continuous form of these 
terms for saddles of any orders. Their coupling constants should be defined from the 
unitarity condition, say in the form of the back-tracking independence of the Wilson 
average (see section 3). 

One possibility to describe the external geometry is to introduce the external 
curvature-dependent terms, like 

' d 2 ^[A(g)X,} 2 (6.2) 

as was proposed by Polyakov. It is not clear whether this simple introduction of the 
extrinsic geometry dependence is sufficient to describe the same effects as saddles do 
for QCD. 

The saddles seem to be relevant factors for the whole problem. In 2 dimensions, 
they define the so-called branch points of the world sheet |7|, |16, [18|, [T?| , like on the 



Riemann surface of the function y = ^fx. This branch point can have a negative 
weight, as seen from the expression for the Wilson average with the loop C encircling 



twice some two dimensional domain of the area A 21 



W(A) = e~ p A {1 - p~ x A) (6.3) 

The last term in the pre-exponent emerges from the entropy of the branch point 
on the minimal area surface covering this contour. The weight of the branch point 
is negative, and for A > (3 the Wilson loop becomes negative. One might wonder 
whether the same phenomenon could happen in the four dimensional QCD in the deep 



confinement regime, when only the minimal area gives the main contribution. For 
the same contour the minimal surface should not look different as in two dimensions, 
with the same branch point. 

The only way to check it is by a Monte-Carlo simulation for the lattice gauge 
theory. A good indication of the relevance of special saddle points on the world sheet 
of the QCD string would be the negative sign of this Wilson average. It would be a 
completely non-perturbative result since in perturbation theory the Wilson average 
is always close to one. Work on this calculation is in progress on the APE in Rome 

Of course, one can easily imagine a much more pessimistic scenario. The saddles 
could be quite big and densely distributed on a typical world sheet for the physically 
relevant values of (3. Such a condensation of saddles might completely cripple the 
surface, and the whole surface description might be irrelevant in this situation. This 
looks to be a probable scenario for big (3, as we have already shown in the one 
plaquette model, where big saddles dominated in the simple sum over "surfaces". 
However, one might hope that these saddles are sufficiently rare on the surface in the 
deep confinement regime. If the worst scenario turns out to be true, it might mean 
that our hopes for a string representation of QCD are hopeless. 
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Figure Captions 

Fig.l. The typical strong coupling diagram for the vector field. 

Fig. 2. Glueing together of two plaquettes after a link variable integration. 

Fig. 3. Formation of saddles, tubes et.c. from four plaquettes after a link variable 
integration. 

Fig. 4 Multiple saddle. 

Fig. 5 Backtracking condition for the wilson loop, following from the unitarity. 



